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We examine several recent lattice-simulation data sets, asking whether they are consistent with 
infrared conformality. We observe, in particular, that for an SU(3) gauge theory with 12 Dirac 
fermions in the fundamental representation, recent simulation data can be described assuming in- 
frared conformality. Lattice simulations include a fermion mass m which is then extrapolated to 
zero, and we note that this data can be fit by a small-m expansion, allowing a controlled extrap- 
olation. We also note that the conformal hypothesis does not work well for two theories that are 
known or expected to be confining and chirally broken, and that it does work well for another theory 
expected to be infrared conformal. 



PACS numbers: 11.30.Qc,11.25.Hf,11.15.Ha,12.60.Nz 
I. INTRODUCTION 

During the past few years, lattice simulations have 
been employed to study the infrared behavior of a va- 
riety of gauge theories that could be relevant to physics 
beyond the standard model. Since much of the numeri- 
cal code was originally developed for QCD, many of the 
simulations have focused on an SU (3) gauge theory with 
varying numbers of massless fermions. 

The lattice simulations of Refs. [TJ [2] considered both 
8 and 12 massless fermions in the fundamental repre- 
sentation, concluding that the former is a confining and 
chirally broken theory like QCD, and indicating that the 
latter is conformal in the infrared, dominated by a fixed 
point. Since then, various authors have examined the 12- 
fermion theory, some agreeing that it is indeed infrared 
conformal [3H2]j but others arguing that it is confining 
and chirally broken [7j [8]. The SU{3) theory with 10 
massless Dirac fermions in the fundamental representa- 
tion has also been studied, with the conclusion that it is 
infrared conformal [5]. 

The recent study of the 12-fermion SU(3) theory by 
Fodor et al [7J is particularly interesting because the sim- 
ulation data set covers a wide range of fermion mass val- 
ues, and finite-volume effects arc relatively small. Their 
analysis leads them to the conclusion that their simu- 
lation data for masses, the pseudoscalar decay constant 
and the chiral condensate are more compatible with con- 
finement and chiral symmetry breaking in the massless 
limit. 

In this paper, we examine the simulation data of Ref. 
[7J noting that it can also be described assuming that 
the massless theory is conformal in the infrared. Lattice 
data for several other theories are also considered. Since 
lattice simulations are carried out by including a fermion 
mass to that is then extrapolated to zero, a question for 
any fit is whether the lattice data can be interpreted in 
terms of a small-m expansion, allowing for a controlled 
extrapolation to zero. We argue that this is the case 
with the conformal hypothesis. Finite- volume effects are 
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considered and shown to be relatively small. 

As a check on this conclusion, we attempt a similar fit 
to an 577(3) theory with 2 fermions in the fundamental 
representation, which is known to be in the broken phase, 
and an SU(3) theory with 6 fermions in the fundamental 
representation, which is strongly believed to be in the 
broken phase. In each case the quality of the fit is poor. 
We also examine the lattice data of Bursa et al [10] for 
an SU(2) gauge theory with 2 fermions in the adjoint 
representation, which is believed to be conformal in the 
infrared PP - TH] . Like the SU(3) theory with 12 fermions, 
the data can be well fit by the conformal hypothesis, with 
a controlled extrapolation to to = 0. 

II. THE CONFORMAL FRAMEWORK 

We first describe the scaling behavior we use to fit 
the lattice data of Ref. 7 . The discussion is similar 
to that in Refs. [13 HE], except that we also include 
nonlcading terms in the scaling behavior. We assume 
that the infrared fixed point g* approximately governs 
the behavior of the theory below some scale A, which, in 
a lattice setting, we take to be the inverse lattice spacing. 
Finite-volume corrections will be described in Sec. IV. 

An explicit fermion mass, m(A) = to is introduced, 
with to <C A. At scales below A, the running mass takes 
the form 

m(/i)=m(A//i)^, (1) 

where 7* > is the mass anomalous dimension evaluated 
at the fixed point g*. At some scale M <C A, the running 
mass satisfies 

to(M) = M. (2) 

At scales below M, the fermions decouple, and the run- 
ning coupling flows away from the fixed point, triggering 
confinement. If the would-be fixed-point coupling g* is 
reasonably strong, the induced confinement scale is of 
order M. We assume this to be the case. 
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The mass of each physical state X is then set by 
the scale M, That is, using Eqs. 1 and 2, Mx — 
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[1/(1+7*)] 



[T7] . where the masses are expressed in 
units of A, and Cx is a dimensionless coefficient not far 
above unity. In addition, there are correction terms, the 
largest of which in a small-m expansion is of order to. 
Keeping only these two terms, we have 
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C x ml 1 ^ 1 ^*'! +D X m. 



(3) 



Since the explicit breaking of chiral symmetry is of order 
M at the induced confinement scale M, there is no ap- 
proximate chiral symmetry to be broken spontaneously. 
Thus this scaling law applies as well to the pseudoscalar 
state. The exponent [1/(1 +7*)] is universal. 

Fodor et al [7J also compute the pseudoscalar decay 
constant F and the chiral condensate ('t/jip) as a function 
of to. Although F plays no special role in the absence of 
spontaneous chiral symmetry breaking, we include it in 
our fit, using an expression similar to that for the masses: 



F = C F m 



[1/(1+7*)] 



Dpm. 



(4) 



The chiral condensate, defined at the cutoff scale A, 
also vanishes as m — > 0. The leading, small-m term is 
purely ultraviolet. This is the "contact term," propor- 
tional to mA 2 , independent of the form of the RG run- 
ning of the coupling and m(p). The second, coming from 
the RG running of (ipi^) from M to A, is proportional to 
M"(3-7 ■'A 7 . Using Eqs. 1 and 2 to express M in terms 
of m and A, we have 



(ipip) — A(jm + Bern 



[(3- 7 *)/(l+7*)] 



(5) 



where now, as in Eqs. 3 and 4, all dimensionful quantities 
are expressed in terms of A, the inverse lattice spacing. 
The coefficients are dimensionless, and to is the lattice 
mass. 

In addition to these terms, we expect a contribution of 
order M 3 , analogous to the leading-order terms in Mx 
and F, arising from the induced confinement scale M. 
And as with Mx and F, there are further corrections, 



one of which is of order 



We therefore take 



(tpip) = Acm + Bern 



[(3-7*)/(l+7*)] 



D c m 3 . 



(6) 



It will turn out that < 7* < 1, so that these four terms 
also provide the basis for a small-m expansion. 



III. FITTING THE LATTICE DATA 
NEGLECTING THE D TERMS 

We fit the lattice data of Ref. [Tj for the masses of 
the scalar, pseudoscalar, vector, axial vector, nucleon, 
and parity partner of the nucleon, for the pseudoscalar 
decay constant, and for the condensate, first setting the 
IP-term coefficients Dx, Dp, and Dc to zero. We then 
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FIG. 1: Log-log plot showing masses of the pseudoscalar 
(P), vector (V), and nucleon (N) states, the pseudoscalar de- 
cay constant (F), and the condensate (C) as a function of 
m, as reported in Ref. [7], along with our conformal fit to 
these quantities, with the D terms set to zero. The universal 
slope of the P, V, N, and F curves provides a good fit to the 
simulation data. 



ask whether the inclusion of the D terms as well as finite- 
volume corrections improves the quality of the fit. 

The simulations of Ref. [7] were performed using a 
tree-level, Symanzik-improved gauge action, with lattice 
gauge coupling j3 = 6/g 2 = 2.2. We assume here that this 
lattice coupling is consistent with the theory being ap- 
proximately described by the infrared-fixed-point value of 
the running coupling throughout the range M < /i < A. 

The simulations were done for fermion masses m = 
0.035, 0.0325, 0.030, 0.0275, 0.025, 0.020, 0.015, 0.010 
(in lattice units), with lattice volume 24 3 x 48 for the 
heaviest 4 masses, with volume 32 3 x 64 for m = 0.025, 
with volume 40 3 x 80 for m = 0.020, and with three 
volumes ranging up to 48 3 x 96 for m = 0.015 and 0.010. 
In the fits reported here, we use the data at the largest 
volume available at each to value [7J. 

In Fig. [T] we show the simulation data for the pseu- 
doscalar (P), vector (V), and nucleon (N) masses, for the 
pseudoscalar decay condensate (F), and for the conden- 
sate (C) as a function of fermion mass to, along with our 
conformal fit to these quantities. The common log-log 
slope for P, V, N, and F, enforced by the universal scal- 
ing exponent 1/(1 + 7*), fits the data points well. The 
slope of the condensate curve is determined dominantly 
by the leading, linear term of Eq. 6. 

To further explore the conformal fit, we report on the 
left side of Table|T]the results of a fit to all the masses, the 
pseudoscalar decay constant, and the chiral condensate, 
with the D terms set to zero and neglecting finite-volume 
effects. For this fit, the anomalous dimension is 7* ps 
0.386 ±0.010 and x 2 /N = 2.508, with N = 53 degrees of 
freedom. 

Our (statistical) error analyses here and below do not 
take into account correlations of the numerical data, 
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Obs. 


D x = 




7* 


0.3858(98) 


0.403(13) 


C P 


4.445(83) 


4.267(85) 


zp 


0.209(64) 


C s 


5.99(14) 


5.75(14) 


zs 


0.63(45) 


C v 


5.26(10) 


5.05(10) 


zv 


0.319(88) 


Ca 


6.68(15) 


6.41(15) 


ZA 


0.50(30) 


Cn 


8.04(17) 


7.70(17) 


ZN 


0.35(18) 


Cjv* 


8.06(17) 


7.73(17) 


ZN* 


0.49(24) 


C F 


0.692(13) 


0.455(39) 


ZF 


0.61(27) 


D F 


— 


0.61(10) 






Ac 


13.898(28) 


13.926(31) 


zc 


-0.036(43) 


B c 


-50.8(5.5) 


-42.2(4.8) 






C c 


94(11) 


79.0(9.5) 






X 2 /dof 


133/53 


42/44 



TABLE I: For the SU(3) theory with 12 fermions in the fun- 
damental representation, best-fit results to the data of Ref. 
[7], for our global conformal fit as described in the text. On 
the left-hand side, all the D terms are set to zero, and finite- 
volume corrections are neglected. On the right-hand side, Df 
[Eq. 4] is included, as well as finite-volume corrections, with 
zx as explained in Sec. IV. The letters S, P, V, A, N, and N* 
correspond respectively to the scalar, the pseudoscalar, the 
vector, the axial vector, the nucleon, and the parity partner 
of the nucleon. F refers to the pseudoscalar decay constant 
and C to the condensate. For each quantity, there are 8 data 
points, one for each m value. The fits do not take into ac- 
count possible correlations between different observables, as 
discussed in the text. 



which would require access to the full simulation data 
set. As a result, our errors may be underestimated and 
X 2 /N may not directly indicate goodness-of-fit. While 
this value of \ 2 /N is somewhat large, the fit reported so 
far is a very simple one. We have not yet included the 
D terms of Eqs. 3, 4, and 6, and we have not taken 
into account possible finite-volume corrections. From 
the fit reported so far, the latter can be anticipated to 
be relatively small. The product ML lies between 1.73 
and 2.23, and finite-volume corrections should bc small 
if M X L cx ML 3> 1. Inspection of the best-fit C x values 
in the left side of Table U indicates that this is the case. 

In the next section, we will include finite- volume effects 
as well as the D terms. The result, reported in the right- 
hand side of Table |TJ is quite encouraging. The finite- 
volume corrections, while not insignificant, are indeed 
relatively small, and the value of 7* changes very little. 
Furthermore, the fit improves, with x 2 /N = 0.944 and 
N = 44. 

Before including the higher-order corrections, we ex- 
amine the consistency of the conformal fit described 
above by performing separate fits to each mass, as well as 
to F and the condensate. We fix a value of 7* in the range 
< 7* < 1, and plot the \ 2 for each fit as a function of 
7* in this range. The result is shown in Fig.[2j along with 
the sum of the individual x 2 's. The internal consistency 
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FIG. 2: For the SU(3) theory with 12 fermions in the fun- 
damental representation, individual contributions to total % 2 
from each channel, for a range of fixed values < 7* < 1. 
Here, all D terms are set to zero. Fits are to the Nf = 12 
data of Ref. [7]. Parity-odd states (P,A,N*) are shown as 
dashed curves in the same color as their parity partners. The 
black curve is the total \ 2 . 



is evident, with the minimum x 2 for each mass occurring 
at a similar value of 7*, and the minimum for the pseu- 
doscalar decay constant at a value only slightly smaller. 
The condensate makes only a small contribution to the 
overall x 2 ■ The minimum of the black curve corresponds 
to the total x 2 of the left fit of Table [i] 



IV. HIGHER-ORDER TERMS AND 
FINITE- VOLUME EFFECTS 

We next comment on the role of the D terms in Eqs. 
3, 4 and 6, inserting one D term at a time and repeating 
the above global fit. At the same time, we consider the 
effect of finite-volume corrections. We incorporate the 
latter by modifying Eqs. 3 and 4 to include first-order 
corrections in an expansion in 1/ML, continuing to in- 
clude the correction linear in the short-distance mass m. 
We thus take M x = C x M [1 + z x /ML] + D x m, with 
M = m\ 1 ^ 1+1 )1, and similarly for F. For the conden- 
sate, we simply replace M = m^^ 1+1 >' in the second 
and third terms of Eq. 6 by M [1 + z c /ML}. While this 
modification is not claimed to be unique, it is physically 
sensible. Furthermore, the finite-volume corrections are 
found to be insignificant for the condensate, so that the 
overall fit does not depend sensitively on the detailed 
form of this correction term. 

For the data set of Ref. [7J, the inclusion of the D 
term in any one of the masses, with or without the finite- 
volume corrections, does not improve the quality of the 
fit. The value of each D x is consistent with 0, with 
errors such that the D x [next-to-leading order (NLO)] 
term is small compared to the C x [leading order (LO)] 
term for the full range of m values. The inclusion of 
the Dq (NNNLO) term in the condensate also does not 
improve the quality of the fit, with Dq consistent with 
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and the errors such that this term is relatively small 
throughout the range of m values. 

The inclusion of the Dp (NLO) term in the pseu- 
doscalar decay constant does improve the global fit, as 
does the inclusion of the finite-volume correction terms 
zx- In the right-hand side of Table |TJ we report the 
results of a fit including Dp and all zx- The resulting 
anomalous dimension is 7* = 0.403(13). The values of 7* 
and the A, B, and C coefficients change very little from 
the left-hand fit. We find Dp = 0.61 ±0.10, showing that 
the NLO term is small compared to the LO term for the 
full range of m values. The values of the zx , zf, and zc 
coefficients show that finite-volume corrections are rel- 
atively small for the entire range of m values, and ex- 
tremely small for the condensate. Most notably, we find 
X 2 /iV = 0.944 with N = 44 degrees of freedom, again 
neglecting possible correlations between the observables. 

In order to make a direct comparison to the broken- 
symmetry fit of Ref. [7j, we have also applied our con- 
formal fit to just the set of four quantities considered in 
their global analysis: the pseudoscalar mass and decay 
constant, the nuclcon mass, and the chiral condensate. 
We find a best-fit value 7* = 0.414 ± 0.016, and with 
other parameters in the same range as reported in Ta- 
ble [l] Most importantly, we find x 2 / N = L1 with N = 
20, to be compared to their broken-symmetry value of 
X 2 /N = 1.22 with N = 26. 

If a global fit including all the D terms is attempted 
for the data set of Ref. [7], the \ 2 dependence on 7* 
becomes very flat, slightly favoring larger values in the 
range < 7* < 1, and leading to poor determinations of 
the D coefficients, with errors comparable in magnitude 
to the central values. We conclude that the current data 
set is not extensive enough to perform a global fit with 
all the parameters of Eqs. 3, 4 and 6. The availability of 
additional simulation data for larger m values would be 
especially helpful in allowing a global fit that constrains 
the D terms. 

For the condensate, the constants Ac, Be and Cc are 
such that the LO (Ac) term strongly dominates through- 
out the mass range. The NLO (Be) term, also arising at 
scales of order A, and the NNLO (Cc) term, arising at 
scales of order M, are of opposite sign and strongly neg- 
atively correlated, with the NLO term dominating the 
NNLO term except for the largest, m = 0.035 point. 

That the fit leads to an NLO coefficient Be of opposite 
sign to the LO term is perhaps surprising, but we know of 
no reason why this correction cannot be of opposite sign. 
It is also possible that this is a consequence of the lim- 
ited amount of simulation data for the condensate, which 
is strongly dominated by the LO, linear term. Further- 
more, since the condensate is so sensitive to physics at 
the scale A, our assumption that the coupling can be 
approximated by its infrared fixed-point value out to A, 
which determines the form of the Be term, should break 
down first here. 

Finally, it is interesting to note that for the masses, the 
Cx coefficients of Table [T] are such that for small m, none 
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FIG. 3: For the SU(3) theory with 2 and 6 fermions in the 
fundamental representation, individual contributions to total 
X 2 from each channel, for a range of fixed values < 7* < 
1. Total number of degrees of freedom is 27 in both cases. 
All D terms are set to zero. The top (bottom) panel shows 
results from fitting to Nf = 2 (Nf — 6) data, obtained from 
the simulations detailed in [181 [19] . Parity-odd states (P, A) 
are shown as dashed curves in the same color as their parity 
partners. 



of the states can decay into a combination of the others. 
Recalling that there is induced confinement in this theory 
at scale M oc ml 1 ^ 1+7 '1, decay into the fundamental 
fermion and gauge-boson constituents is forbidden. Since 
there is no reason for other states such as the ++ to be 
lighter, it appears that each of these states is stable for 
arbitrarily small m, but with an induced confinement 
radius diverging as m — > 0. 



V. TESTING INFRARED CONFORM ALITY ON 
OTHER GAUGE THEORIES 

As a check on our conclusion that the simulation data 
of Ref. [7| are consistent with infrared conformality, we 
have tested our conformal fit on two theories for which 
it should not work well and one for which it should. The 
former are an SU (3) gauge theory with 2 fermions in the 
fundamental representation, known to be in the broken 
phase, and an SU(3) gauge theory with 6 fermions in 
the fundamental representation, strongly believed to be 
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in the broken phase. Here we fit the simulation data 
of Refs. [HI [19]. In each case the quality of the fit is 
indeed poor, as shown in Fig. [3J plotting the x 2 for each 
individual fit, as well as the overall x 2 > as function of 7*. 
No clear minimum in x 2 appears for any channel except 
the Nf = 6 pseudoscalar mass, where it is at 7* close to 
1. For Nf — 2, a minimum appears at 7* ~ 1. With 
chiral symmetry breaking, Mp ~ m 1 / 2 in lowest order 
chiral perturbation theory, corresponding effectively to 
7* = 1. 

It is worth noting that for the SU(3) gauge theory with 
2 or 6 fermions in the fundamental representation, the 
poorness of the conformal fit should not be due to finite- 
volume effects. Within the conformal hypothesis, as we 
noted in the case of the 12-fermion theory, a measure of 
finite-volume effects is given by ML = m 1 ^ 1+7 ^L. Here, 
the value of 7* emerging from the poor fit is of order 
unity, and L = 32, so ML > 2.3 for the entire range of 
m values. Each of the associated masses is larger than 
M, so finite-volume effects should be relatively small. 

We also note that for the SU(3) gauge theory with 
2 fermions in the fundamental representation, a fit us- 
ing chiral perturbation theory for a confining and chiral- 
breaking theory does work well [20] . For the SU (3) gauge 
theory with 6 fermions in the fundamental representa- 
tion, a smaller set of fermion masses will be required to 
apply chiral perturbation theory |20) . But there is strong 
evidence from lattice simulations of the running coupling 
that this theory is in the broken phase [TJ [5] • 

A theory for which a conformal fit should work well 
is an SU(2) gauge theory with 2 fermions in the adjoint 
representation, widely believed to be conformal in the in- 
frared [TTHT41 . We have fit the simulation data of Bursa 
et al [TU] for the pseudoscalar and vector masses and 
pseudoscalar decay constant, assuming as before that M 
is the induced confinement scale up to a coefficient of or- 
der unity. The data are used only in the range m < 0.2, 
in order to ensure that our formulas based on a small-m 
expansion can be applied. Figure [4] shows x 2 versus 7* 
for each channel, as well as the overall x 2 j based on con- 
formal fits as in Eqs. 3 and 4. We find a clear minimum 
in x 2 at the best-fit value of 7* = 0. 17 ± 0.05. This value 
is roughly consistent with previous determinations of 7* 
[TTl [TJ]. Again, our analysis does not include the full 
data covariance matrix, so our error may be underesti- 
mated and x 2 /-V may not indicate goodness-of-fit. The 
relatively large contribution of the decay constant to the 
overall x 2 ma y be due to underestimation of statistical 
errors, as the data as shown in Ref. [TU] are difficult 
to describe with any smooth function of m. Addition 
of finite-volume corrections as described in Sec. IV sup- 
ports this conclusion, failing to improve the quality of 
the decay constant fit (but improving fits to the masses.) 
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FIG. 4: For the SU (2) theory with two fermions in the adjoint 
representation, individual contributions to total \ 2 from each 
channel, for a range of fixed values < 7* < 1. Total number 
of degrees of freedom is 14 (6 data points per channel). All 
D terms are fixed to zero. Fits are to the SU(2) two-flavor 
adjoint data of Bursa et al [10], with the restriction m < 0.2. 



VI. CONCLUDING COMMENTS 

We have argued that the simulation data of Ref. [7] 
are consistent with the hypothesis that an 5/7(3) gauge 
theory with 12 massless fermions in the fundamental rep- 
resentation is conformal in the infrared. This conclusion 
is based on a simple fit to the data, in particular assum- 
ing that the gauge coupling can be approximated by its 
infrared-fixed-point value g* out to the UV cutoff A (the 
inverse lattice spacing). The mass anomalous dimension 
7 is then set to its fixed-point value 7*. 

A global fit including finite-volume and higher-order 
corrections yields x 2 /N — 0.944 with N — 44. We have 
argued that the fit describes a small-m expansion cov- 
ering the range of fermion masses used in the simula- 
tions, allowing a controlled extrapolation to m = 0. It 
leads to a mass anomalous dimension of 7* = 0.403(13). 
To compare directly with the broken-symmetry fit of 
Ref. [7], we also fit to a subset of four channels, finding 
X 2 /A = 1.10 with N = 20 as compared to the reported 
broken-symmetry value of X 2 /N = 1.22 with N = 26. 

Although not described in detail here, we have also 
used the infrared-conformal hypothesis to fit the simula- 
tion data of Ref. [7] for the static quark potential. Since 
confinement is induced at scale M, an effective string 
tension a cx M 2 oc m\ 2 ^ 1+1 " can be determined from 
the data assuming that string breaking has not yet set in. 
The fit works well, with a value of 7* in good agreement 
with the other fits and an acceptable x 2 - 

We stress that we have not argued conclusively that the 
simulation data of Ref. [7] demonstrates that the SU(3) 
theory with 12 massless fermions is infrared conformal. 
The simulation data can be described with similar fit 
quality by the chirally broken functional forms used in 
7J, with a slope and intercept. But the large value of the 
slope term compared to the intercept, for the existing 
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range of m values, does not provide the basis for a small- 
m expansion in the spirit of chiral perturbation theory 
with a controlled extrapolation to m = 0. Further simu- 
lations at additional m values will help to distinguish the 
two scenarios. 

As a check on the validity of our infrared-conformal 
fit to the SU(3) theory with 12 fcrmions, we have ob- 
served that it does not work well for an SU(3) gauge 
theory with 2 fermions in the fundamental representa- 
tion, which is known to be in the broken phase, or an 
SU(3) gauge theory with 6 fermions in the fundamental 
representation, which is strongly believed to be in the 
broken phase. On the other hand, it does work well for 
an SU (2) gauge theory with 2 fermions in the adjoint 
representation, believed to be in the infrared-conformal 
phase. We have fit the pseudoscalar and vector masses 
and decay constants, with a 7* consistent with other ref- 
erences. Here, too, we have fit the string tension induced 
at finite m with an anomalous dimension in agreement 
with the other fits. 

It will next be important to address several of the sim- 
plifying assumptions made in these fits. Finite-volume 
effects should be examined further, as well as corrections 



due to the running of the coupling and mass anomalous 
dimension at higher mass scales. Also, a possible hierar- 
chy between M and the induced confinement scale should 
be considered. 
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